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Remnants of massive metal-poor stars: viable engines for ultra- 
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Massive metal-poor stars might end their life by directly collapsing into massive (« 25 — 80 Mq) black holes (BHs). We 
derive the number of massive BHs (Nbh) that are expected to form per galaxy via this mechanism. We select a sample 
of 66 galaxies with X-ray coverage, measurements of the star formation rate (SFR) and of the metallicity. We find that 
Nbh correlates with the number of observed ultra-luminous X-ray sources (ULXs) per galaxy (Nulx) in this sample. We 
discuss the dependence of Nulx and of Nbh on the SFR and on the metallicity. 
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1 Introduction 



Most ultra-luminous X-ray sources (ULXs) are located in 
galaxies with a high star formation rate (SFR, e.g. Irwin, 
Bregman & Athey 2004). The ULXs match the correlation 
between X-ray luminosity and SFR reported by various stud- 
ies (Grimm, Gilfanov & Sunyaev 2003; Ranalli, Comastri & 
Setti 2003; Gilfanov, Grimm & Sunyaev 2004a,b,c; Kaaret 
& Alonso-Herrero 2008; Mineo & Gilfanov 2010). Further- 
more, the same studies indicate that the luminosity function 
of ULXs is the direct extension of the function for high- 
mass X-ray binaries (HMXBs). Recent papers suggest a cor- 
relation between ULXs and low-metallicity environments, 
and propose that this may be connected with the influence of 
metallicity on the evolution of massive stars (Pakull & Miri- 
oni 2002; Zampieri et al. 2004; Soria et al. 2005; Swartz, 
Soria & Tennant 2008). This scenario has been explored in 
detail by Mapelli, Colpi & Zampieri (2009, hereafter M09), 
by Zampieri & Roberts (2009) and by Mapelli et al. (2010, 
hereafter M10), highlighting that a large fraction of ULXs 
may actually host massive (~ 30 — 80 Mq) stellar black 
holes (BHs) formed in a low-metallicity environment. In 
fact, low-metallicity (Z ^ 0.4 Z@) massive stars lose only 
a small fraction of their mass due to stellar winds (Maeder 
1992, hereafter M92; Heger & Woosley 2002, hereafter HW02: 
Heger et al. 2003, hereafter H03; Belczynski et al. 2010, 
hereafter B10) and can directly collapse (Fryer 1999; B10) 
into massive BHs (25 M Q < ttibh < 80 M Q ). These mas- 
sive BHs can power most of the known ULXs without re- 
quiring super-Eddington accretion or anisotropic emission. 
Furthermore, their formation mechanism can explain the 



correlation between ULXs and SFR, and the fact that ULXs 
are preferentially found in low-metallicity regions. 

2 Sample of galaxies 

In this proceeding, we consider a sample of 66 galaxies. 
All of them have X-ray coverage, at least one measurement 
of the star formation rate (SFR) and of the metallicity (Z). 
64 galaxies are taken from the sample listed in Table 1 of 
M10. The remaining two are I Zw 18 and the interacting 
pair SBS 0335052/SBS 0335052wfl These two objects are 
extremely metal-poor galaxies (XMDs, Moiseev, Pustilnik 
& Kniazev 2010, and references therein) and are impor- 
tant, because they are the only galaxies with Z < 0.03 Zq 
and with X-ray observations. Their properties and the corre- 
sponding references are listed in Table 1 . For details about 
the data and the properties of the other 64 galaxies, see M10. 

For all the galaxies in the sample, we derive a fidu- 
cial value for the SFR (when there is more than one mea- 
surement, we take, in general, the average value), for the 
metallicity (we adopt an uniform calibration, see Pilyugin 
& Thuan 2005; when a metallicity gradient is available, we 
take the value of Z at 0.7 Holmberg radii, see M10) and we 
estimate the number of ULXs Ntjlx after subtracting the 
background contamination (see M10 for details). 

3 Observational results 

The data collected from the literature were analyzed follow- 
ing the same procedure as described in M10. In particular, 
we adopt the \ 2 analysis (although such method might not 
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1 In this proceeding, we consider the interacting pairs as a unique object, 
for consistency with M10. 
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Table 1 Properties of the two XMDs in the sample. 



Galaxy 


SFR[M vr" 1 ]* 


Z [Z e ] b 


N ULX C 


IZw 18 


0.07 


0.02 


1 


SBS 0335-052 


1.1 


0.025 


2.83 



a SFR from Wu et al. (2007) for I Zw 18 and from Thuan, Izotov & 
Lipovetsky (1997), Pustilnik, Pramskij & Kniazev (2004) and Johnson, 
Hunt, Reines (2009) for SBS 0335-052. b Metallicity from Thuan et al. 
(2004) for both galaxies. c Nulx is the number of ULXs per galaxy after 
subtracting the background contamination. Thuan et al. (2004) find 1 
ULX in I Zw 18 and 3 ULXs in SBS 0335-052/SBS 0335-052W. The 
estimated contamination is 0.00 in I Zw 18 and 0.17 (upper limit) in 
SBS 0335-052/SBS 0335-052W. 
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Fig. 1 Upper panel: Nulx versus the SFR. Filled black 
circles: galaxies with metallicity < 0.2 Zq\ open circles 
(red on the web): galaxies with metallicity > 0.2 Zq. Solid 
line: power-law fit for the entire sample; dashed line (red 
on the web): power-law fit obtained assuming that the index 
of the power law is equal to 1. Central panel: Nulx versus 
Z. Filled black circles: entire sample; solid line: power-law 
fit. Lower panel: Nulx/SFR versus Z. Filled black circles: 
entire sample. Solid line: power-law fit. In all the panels: the 
error bars on both the x— and the y— axis are 1 a errors. 



be completely suitable for small samples, see M10), per- 
form the power-law fits and calculate the correlation co- 
efficients (see Table 2). From the upper panel of Fig. Q]it 
appears that there is a strong correlation between Nulx 
and the SFR in our sample. Such correlation is consistent 
with a linear relation (see the best-fitting values reported 
in Table 2), in agreement with previous studies (see e.g. 
Grimm, Gilfanov & Sunyaev 2003). Instead, no significant 
correlation appears between Nulx an d Z (central panel of 
Fig. [T]). However, a marginally significant correlation ex- 



ists between the number of ULXs normalized to the SFR 
(Nulx/SFR) and the metallicity. This suggests that the metal- 
licity affects the formation of ULXs, but its contribution is 
less important than that of the SFR. The idea that metal- 
licity plays a role in the origin of ULXs is consistent with 
previous observations (see e.g. Swartz et al. 2008 and ref- 
erences therein) and with some recent theoretical models 
(M09; Zampieri & Roberts 2009; Linden et al. 2010; M10). 
Finally, the presence of the two XMDs in our sample does 
not significantly change the best-fitting values (Table 2) with 
respect to those derived in M10. 

4 Comparison of the data with the 
theoretical model 

In this Section, we analyze the observational data collected 
from the literature on the light of the theoretical model re- 
cently proposed by M09 and M10. First, we briefly summa- 
rize such model. 

4.1 Theoretical model 

According to numerical calculations (Fryer 1999; HW02; 
H03), a star that, at the end of its life, has a final mass 
m&n > 40 M is expected to directly collapse into a BH. 
In this case, the mass of the remnant BH is likely more 
than half of the final mass of the progenitor star, as rela- 
tively small mass ejection is expected in the direct collapse. 
Thus, stars that at the end of their lives have mfi n > 40 M© 
are likely to produce massive BHs (B10). The final masses 
of the stars strongly depend on their metallicity. Massive 
stars with metallicity close to solar cannot have final masses 
larger than mfj n ~ 10—15 M , even if their initial mass was 
very large, as they are expected to lose a lot of mass due to 
stellar winds (H03). Instead, massive stars with lower metal- 
licity are less affected by stellar winds, and retain a larger 
fraction of their initial mass. If its metallicity is sufficiently 
low, a star can have a final mass > 40 M and can 
directly collapse into a massive BH with a mass 25 M Q < 
m B H < 80 M Q (HW02; H03; B10). The recent model by 
B10 accounts for the fact that stars with ma n > 40 M Q can 
directly collapse into BHs. For this reason, in this model 
BHs with mass as large as 80 M Q are allowed to form. 

On the basis of this scenario, we can derive the expected 
number of massive BHs per galaxy (Nbh) as a function of 
the star formation rate (SFR) and of the metallicity Z (see 
M09, M10): 

/>m max 

Nbh (SFR, Z) = A(SFR) / m~ a dm, (1) 

Jm proB (Z) 

where m max is the maximum stellar mass (we assume m max = 
120 M Q ) and a is the slope of the initial mass function 
(IMF). Here, we adopt the Kroupa IMF, for which a = 1.3 
if m < 0.5 Mq and a = 2.3 for larger masses (Kroupa 
2001). m pr0 g(Z) is the minimum initial stellar mass (i.e. the 
mass at zero-age main sequence) for which a star is the pro- 
genitor of a massive BH. As we discussed above, m plog (Z) 
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Fig. 2 Upper panel: Nbh (derived using the model from 
BIO) versus the SFR. Central panel: Nbh versus Z. Lower 
panel: Nbh/SFR versus Z. In all the panels, filled black 
circles: entire sample; solid line: power-law fit for the entire 
sample; the error bars on both the x— and the y— axis are 
1 a errors. 



strongly depends on the metallicity. In our calculations, we 
assume m prog (Z) to be the initial stellar mass for which 
the mass of the remnant is mBH ~ 25 M©, according to the 
model by B 1(0. 

Finally, A(SFR), the normalization constant in equa- 
tion (Q]), can be estimated as 

A(SFR) = rmn f R i° - , (2) 



J ri 



~ a dr 



where m m i n is the minimum stellar mass (we assume TO m ; n = 
0.08 Mq), SFR is the current star formation rate and t co is 
the characteristic lifetime of a possible companion of the 
massive BH. In fact, we are not interested in all the massive 
BHs, but only in those that could acquire massive stellar 
companions and power observable ULXs. In this proceed- 
ing, we adopt a constant value i co = 10 7 yr, which is the 
lifetime of a ~ 15 M Q star. 

4.2 Results 

Fig. |2] shows the behaviour of the theoretical model when 
applied to the observed SFR and metallicity. As assumed 
in the model, Nbh scales linearly with the SFR. The cen- 
tral panel of Fig. [2] shows that there is no significant cor- 
relation between Nbh an d Z, although we imposed, in the 




MK) also consider an alternative model by Portinari, Chiosi & Bressan 



1998. 



Fig. 3 Number of observed ULXs per galaxy Nulx ver- 
sus the number of expected massive BHs per galaxy Nbh, 
derived using the model from B10. The solid line is the 
power-law fit for the entire sample. The dashed line (red 
on the web) is the power-law fit obtained assuming that the 
index of the power law is = 1. The error bars on both the 
x— and the y— axis are 1 a errors. 



model, that Nbh does depend on Z. This absence of correla- 
tion agrees with what we found for Nulx versus Z (central 
panel of Fig. Q]). Finally, only when the effect of the SFR 
is subtracted (by normalizing Nbh to the SFR, lower panel 
of Fig. [3}, it is possible to see the dependence of Nbh on 
the metallicity. The behaviour of Nbh/SFR versus Z in the 
model is consistent with that of Nulx /SFR versus Z in the 
data. 

Fig. [3] and Table 2 indicate that there is a correlation 
between Nbh and Nulx- This correlation is slightly more 
significant than that between Nulx and SFR, when consid- 
ering both the x 2 analysis and the correlation coefficient. 

Recently, Linden et al. (2010) proposed a different model 
to explain the connection between low-metallicity environ- 
ments and ULXs. They indicate that the number, the life- 
time and (less significantly) the luminosity of HMXBs are 
enhanced by low metallicity. Linden et al. (2010) also point 
out a possible problem of MIO's model: massive BHs born 
via direct collapse likely do not receive a natal kick and 
this fact excludes, in the model by Linden et al. (2010), the 
possibility of forming a HMXB via Roche lobe overflow 
(RLO). 

On the other hand, Linden et al. (2010) always require 
super-Eddington emission, to explain the ULXs. Further- 
more, the process of the direct collapse and the physics of 
the binaries which include massive BHs born from it are still 
far from being understood. For example, natal kicks might 
still be present, due to asymmetries induced by sterile neu- 
trinos (e.g. Kusenko 2006). In alternative, kicks might occur 
for different reasons, e.g. due to three-body encounters in 
the parent cluster (Mapelli et al., in preparation). Such sce- 
nario might also explain why ULXs are often found to be 
displaced with respect to star-forming regions (e.g. Swartz, 
Tennant & Soria 2009). 

5 Conclusions 

In this proceeding, we considered a sample of 66 galaxies. 
All of them have X-ray coverage, at least one measurement 
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Table 2 Parameters of the power-law fits and \ 2 



X 


y 


Index a 


Normalization 


x 2h 


r c 


Nbh 
Nbh 


Nulx 
Nulx 


0.82^:2 
1.00 


-2 63 +a48 

Q QJ- + 0.06 

-J.JJ_0.06 


12.0 
13.0 


0.92 
0.92 


SFR 
SFR 


Nulx 
Nulx 


1.00 


o ifi+ 09 

u - iD -0.12 

o nq +0 ' 06 


19.3 
19.8 


0.88 
0.88 


Z 
Z 


Nulx 

Nulx/SFR 


-0 16 +a2S 

O-lO_0.28 
-0 <jS+ a21 


12+ ' 20 

u - lz -0.20 

-0 37+ - 16 

u --"-0.16 


86.5 
11.0 


-0.15 
-0.38 



The SFR and the Z used by the fitting procedure are in units of Mq yr 
and of Zq, respectively. a When the index is equal to 1.00 or to 0.00, 
without error, it means that it has been fixed to such value. b x 2 is the 
non-reduced \ 2 . The number of degrees of freedom (dof) is 65 when the 
index has been fixed, 64 in the other cases. c r is the Pearson correlation 
coefficient. 



of the SFR and of Z. This sample includes two XMDs, 
which have extremely low metallicity and host a relatively 
high number of ULXs, when compared to their SFR. 

We find that there is a strong correlation between the 
number of ULXs per galaxy (Ntjlx) and the SFR. This is 
consistent with previous studies (e.g. Grimm, Gilfanov & 
Sunyaev 2003). We also find a marginally significant anti- 
correlation between Ntjlx/SFR and the metallicity. This 
might indicate that the metallicity is the missing ingredi- 
ent, to understand the formation of ULXs, although the error 
bars are still very large and the sample of galaxies is quite 
small. 

Recently, M09 and M10 suggested that ULXs might 
be connected with massive BHs (25 — 80 Mq) formed by 
the direct collapse of massive metal-poor stars (Fryer 1999; 
B10). We derive the number of BHs per galaxy (Nbh) pre- 
dicted by M10 and compare it with the observed Ntjlx in 
our sample. We find a strong correlation between Nbh and 
Nulx- 

We note that the model by B10 derives the mass of the 
remnant for single stars only, without considering stars in 
binaries. Stars in close binaries likely have a different mass- 
loss history. Accounting for binary progenitors might strengthenrh u 
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the dependence of the BH mass on metallicity. Therefore, it 
will be necessary to account for binary evolution, to refine 
the model of massive BH formation. 

Furthermore, the physics of the direct collapse and the 
properties of massive BHs born from it are only poorly known. 
These need to be investigated, in order to understand the 
process of mass transfer (and of X-ray emission) in binaries 
including massive BHs. 

Finally, we need more observational data, especially mea- 
surements of the metallicity, to strengthen our conclusions. 
XMDs are particularly interesting, because of their pecu- 
liarly low metallicity. 
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